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Abstract 

Assuming that the phase transition between the Wigner solid and the Laugh- 
lin liquid is first-order, we compare ground-state energies to find features of 
the phase diagram at fixed v. Rather than use the Coulomb interaction, we 
calculate the effective interaction in a square quantum well, and fit the results 
to a model interaction with length parameter A roughly proportional to the 
width of the well. We find a transition to the Wigner solid phase at high 
density in very wide wells, driven by the softening of the interaction at short 
distances, as well as the more well-known transition to the Wigner solid at 
low density, driven by Landau-level mixing. 
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In the past few years, experiments with 2- dimensional electron systems in a perpendicular 
magnetic field have shown a long-expected behavior as the strength of the field increases. 
At comparatively low fields, the now-familiar fractional quantum Hall effect appears, as the 
longitudinal resistivity p xx falls exponentially to zero as the temperature T — > at fractional 
filling factors v. As the field increases, however, an insulating state appears, with p xx — > oo 
when T — > 0. HHI! Optical experiments {||| show a new spectral line developing at the same 
fields. This behavior is commonly thought to signal the presence of a Wigner solid, predicted 
many years ago, and pieces of an experiment ally- derived phase diagram between the Wigner 
solid and the fractional quantum Hall (FQHE) or Laughlin liquid can be sketched out. 

There are several experimental parameters that affect the boundary between Wigner solid 
and FQHE liquid. 0-|9[ The most important is the magnetic filling factor v = 2irn£ 2 , where 
n is the electron density and the magnetic length i = (hc/eB) 1 ^ 2 . The FQHE appears only 
at certain fractional filling factors v = p/q, where p, q are integers and q is odd. The Wigner 
solid exhibits no such detailed dependence on u, but becomes gradually more favorable 
as the particles are localized with decreasing v. Also important is the electron density n, 
parameterized by the ion-disk radius r s = {jin)^ 1 / 2 . (Here we use atomic units, where length 
is measured in units of the Bohr radius ag — Ti 2 /me 2 and energy in units of e 2 /a^.) As 
the Landau-level separation is Tiu c = 2/isr 2 , the energy cost of localizing the particles by 
Landau-level mixing falls rapidly with increasing r s at fixed u, until at some critical r s the 
Wigner solid becomes more favorable than the FQHE liquid and the system freezes. 

Finally, the experiments are all done in real systems, which must be considered quasi- 
two-dimensional, with some finite thickness L characterizing the width of the electron layer. 
On average a pair of electrons is separated in the z-direction by A < L, so their effective 
interaction at distances r C A in the xy plane becomes much softer than the Coulomb 
interaction, while at large distances r ^> A, the interaction is essentially Coulombic. This 
preserves the long-range character of the interaction while weakening the short-range part. 
Because the formation of the Wigner solid is driven by the long-range part of the interaction, 
while the FQHE liquid derives its energy advantage from the short-range part, the quasi- 
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two-dimensional character of a real experimental system might be expected to favor the 



solid somewhat. There have been several recent experiments and theoretical studies pO|jTT 
of finite thickness effects on the incompressible FQHE state and, in general, it is now well- 
established that finite layer thickness tends to destroy incompressibility by reducing the 
short-range part of the Coulomb interaction. To the best of our knowledge, however, the 
enhancement of the Wigner solid phase by the finite thickness effect has not been examined 
theoretically. 

In this paper, we study the quantitative effects of the three length scales £, r s , and 
particularly A on the Laughlin liquid - Wigner solid phase diagram. We compare ground- 
state energies and excited states of variational wavefunctions for the Wigner solid and the 
FQHE liquid as both r s and A are allowed to vary. We find that at small A, such as that 
found in GaAs heterojunctions, the zero-temperature phase transition from liquid to solid 
at large r s does not vary much from the Coulomb case. || At large A we find that the FQHE 
liquid gives way to the Wigner solid at low density, as expected, but most unexpectedly, we 
find that the Wigner solid phase dominates the FQHE liquid at high density as well. Only 
in an intermediate range of r s does the liquid have lower energy than the solid when A is 
large. The interplay among the length scales £, r s , and A can, in principle, therefore lead to 
a reentrant Wigner solid transition at large A as r s is varied. 

The electrons in a quantum well are confined to a small number of subbands, usually 
just one, and then can be thought of as extended rod-like charges in the z-direction which 
are allowed to move only in the xy plane. [0 This approximation has been used in the past 
to study the weakening and eventual collapse of the FQHE state in a wide quantum well. 
lT| The interaction between these model charges is given by 



J J [r 2 + [z\ - z 2 y\ 
Here ((z) is the envelope wavefunction describing quantization in the z-direction and r is 
the separation between electrons in the xy plane. 

The confining potential in the z-direction enters into this equation only through the 



envelope wavefunction ((z). This wavefunction should be obtained in a self-consistent pro- 
cedure which takes into account the interaction of the electrons in the xy plane. At low 
electron density, in quantum wells where the subband splitting is large, ((z) is simply the 
z-component of the single-particle wavefunction. When electron density becomes higher, or 
the well is made wider bringing the subbands closer together, ((z) will be modified. 

In a square well, however, A and r s are (roughly speaking) independent parameters, as 
the width of the well L is to lowest order independent of the density of the electrons in the 
xy plane. It is only when subband mixing begins to become important that the electron 
density begins to affect the envelope wavefunctions ({z). We expect that when only the 
lowest subband is occupied the single-particle wavefunctions will be adequate for ((z). The 
usual cosine solution for an infinite well is shown in the inset to Figure [I], as we U as a 
gaussian wavefunction 
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fitted to the cosine solution. Here 7 = 0.277L gives the gaussian wavefunction shown in 
Figure |1|. The effective potential for the two wavefunctions is almost identical. The gaussian 
wavefunction gives 



(3) 

where K$ is a modified Bessel function. We would like to use the simpler model interaction 



of Zhang and Das Sarma, [13| 



Mr) - 7WW , (4) 

so we choose the parameter A to fit V^(r) best when r is large. A least-squares fit, shown 
in Figure 0, yields X/L — 0.2. 

Because X/L is small for the square well, we need a wide well if we are to investigate a 
system with reasonably large A. For example, a well with A = 1 in an electron system in GaAs 
is approximately 500Awide. In these wide wells, subband mixing can become important, 



and the envelope wavefunctions will tend to spread out toward the edges of the well as the 
electrons reduce their potential energy. We have used the self-consistent approach taken in 
|I(J to estimate the density distribution n(z), and we find that X/L will vary anywhere from 
0.15 to 0.30 in the presence of subband mixing. [[H]] If we choose the lowest subband value 
A = 0.2L for the square well, Vo(r) is nearly identical to V^(r) for r > 0.2L, and differs 
significantly from V^(r) only for r < 0.1L. The pair correlation function for both liquid 
and solid is small below about r = r s , so we believe Vo(r) to be a good description of (r) 
for A < 2. 

In extremely wide quantum wells, n(z) will be modified further as it peaks near the well 
edges. The well then begins to resemble a highly-coupled double-layer system. Electrons 
in the well may lower their potential energy by localizing near the edges of the well, at 
the cost of some kinetic energy. In this more complicated situation, the approximation (|]) 
begins to break down. At high density, the interaction between electrons that are adjacent 
in the x|/-plane becomes small, as they become separated in the ^-direction. The interaction 
between electrons localized on the same side of the well becomes stronger, however, and the 
net effect is a small potential energy savings. We should note, however, that the magnetic 
field will tend to suppress subband mixing, as at v = 1/3 the filling in the lowest subband 
will be at most 1/3. 

In this paper, we determine the zero-temperature phase boundary between Wigner solid 
and Laughlin liquid by comparing ground-state energies of variational wavefunctions for 
the liquid and the solid. Because we are simply comparing energies, we are assuming that 
the phase transition is first-order, and we neglect the possible presence of any other states 
in the vicinity. At high electron density we regard the approximation (|]) as a qualitative 
guide only, and do not attempt to predict a critical r s and A quantitatively for the low r s 
transition. 

A variational wavefunction for the liquid which interpolates in some sense between a 
wavefunction with the lowest possible kinetic energy, the Laughlin wavefunction, and a 
wavefunction with the lowest possible potential energy, in which the electrons are completely 
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localized, might be expected to be a good variational choice. To that end, we have chosen a 
variational wavefunction that consists of the Laughlin wavefunction ip m (m = multiplied 
by a Jastrow factor ]li<j e ~ a ^ where is the distance between the i th and j th particles 
and a is the variational parameter. When a = 0, we recover the Laughlin wavefunction, 
and when a^O, the wavefunction is no longer analytic and higher Landau levels are mixed 
in. The Jastrow factor introduces more correlations into the wavefunction, lowering the 
potential energy, while introducing a kinetic energy cost. 

Details of the calculation are given in Reff| so we will only review it briefly here. (Note 
that H uses units of energy e 2 /2ae, while this paper uses atomic units e 2 /ae.) We use the 
spherical geometry, in which our wavefunction becomes 

C = II Mi - u i v i) m ex P 1— ; r - 11/2 ' ( 5 ) 



where u% = e - * 9 ^ 2 cos(0j/2), V{ = e"^ ;//2 sin(#j/2) are convenient spinor coordinates, and the 
distance between particles i and j is taken as the chord distance ry = 2R\uiVj — UjVi\. 
Evaluating the energy of this wavefunction by Monte Carlo and minimizing at a fixed r s 
and A gives the results shown by the dashed lines in Figures and [3]. 

In order to find the liquid-solid phase boundary we need a rather accurate evaluation of 
the solid. Lam and Girvin evaluated the energy of a correlated Wigner solid wavefunction 

* = ex p f iE'&%& ) iiM*), (6) 

V i,3 I i 



where Ci — z i~ Ri, F>ij = B{Ri — Rj)i an d 
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Here Z{ = Xi + iyi is the i th particle position and Ri = Xi + iY^ is the i th lattice site. ^ is 
the harmonic crystal wavefunction restricted to the lowest Landau level, and the variational 
parameters Bij are calculated by using the values derived from the harmonic crystal. How- 
ever, in order to make a reasonable comparison of solid and liquid wavef unctions, we need, 
as discussed above, a wavefunction which includes Landau-level mixing. 



A study of the ground-state energy of the Wigner crystal including Landau-level 
mixing, has recently been completed, using the Coulomb interaction. We have extended 
their work to use the modified potential (Q). This calculation is similar to Lam and Girvin's, 
except that two more variational parameters a and (3 are added to the wavefunction to put 
in more correlations at the expense of some Landau- level mixing. First, the gaussians in the 
single-particle wavefunctions (0) were "squeezed" to move the electrons farther away from 
each other, by making the replacement 

exp [~\ Zi - -> exp (-01* - Ri\ 2 ) . (8) 

Varying the parameter f3 away from 1/4 introduces Landau- level mixing into the wave- 
function because the single-particle wavefunctions <f%.(zi) are no longer eigenstates of the 
single-particle Hamiltonian. An additional Jastrow factor is then introduced, and the final 
wavefunction, with two variational parameters a and /3, is then 

* = exp [£' [i&flyfc - ~t*(|* - ^1)1 ) n<(*), (9) 



where 
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and F is a constant chosen to optimize the pseudo-potential at small r. Zhu and Louie 
varied B^ as well as /3, but found that varying B^ had very little effect on the energy, as 
Lam and Girvin had suggested. The quantum Monte Carlo calculations were done with 
modified periodic boundary conditions, |7| which require the addition of a phase factor to 
the single-particle wavefunctions (|7|). The results are shown as the solid lines in Figure |2| 
and the triangular data points in Figure 

The r s = 2 curves in Figure |2] are very nearly the lowest Landau-level energies, since hu c 
is large at r s = 2. In fact, if each of the curves in Figure were plotted as a function of 
X/r s and there were no Landau- level mixing, they would lie on top of each other. Only the 
increased Landau- level mixing at r s = 10 and 20 reduces the energy there somewhat. The 



parameter X/r s can be thought of as the ratio of the average separation in the ^-direction 
A and the average separation in the xy plane r s (actually ~ 2r s ) of two nearby electrons. 
The lowest Landau level energies cross at A = 6 at r s = 10 and A = 12 at r s = 20, but the 
variational energies here predict a freezing transition at A ~ 4 at r s = 10 and possibly A ~ 6 
- 7 at r s = 20. These results illustrate the fact that Landau-level mixing does in general 
tend to favor the solid somewhat. 

The variational wavefunctions use the fact that localizing the electrons by mixing in 
higher Landau levels will keep nearest-neighbor electrons farther apart, and since the 
Coulomb potential rises rapidly at small r the energy savings can be significant. As the 
interaction softens, however, the energy savings becomes small, and as a result the liquid 
and solid energies at r s = 2 are nearly identical to the lowest Landau level energies. Ortiz, 
Ceperley, and Martin [|15[ have recently used their fixed-phase quantum Monte Carlo method 



to improve the liquid energy shown here somewhat. The difference in energy is significant 
at r s = 10 and 20, but at r s = 2 the change in energy for the Coulomb interaction is only 
~ 0.006, a shift on the order of 10%. 

The rapid descent of the solid energy with A leads to a freezing transition at about A = 1.2 
at r s = 2, brought on by the softening of the short-range part of the interaction. Keeping the 
well width constant and lowering the density gives a smaller X/r s , and the system melts once 
again until the increased Landau-level mixing at high r s causes the system to freeze again. 
Figure ^| shows the energies of solid and liquid when A is fixed and r s is allowed to vary. At 
A = 0.5 we find no significant lowering of the solid energy at low r s , but at A = 1 the solid 
energy at r s = 2 is nearly as low as the liquid energy, and at A = 2 and 3 the solid energy 
at r s = 2 is much lower than the liquid energy. In all cases the liquid becomes favorable at 
lower density until Landau-level mixing again causes the system to freeze. The transition at 



high density has been seen in the experiment of Suen et al, |TIJ where an insulating phase 
was observed in an 800Awell at v — 1/3, r s = 1.7, which became a well-defined FQHE state 
when the density was lowered to r s = 2.2. 

We can contrast the wide well (i.e. large A) situation, which we have argued in this 
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paper may be favorable to the solid phase both at large and smaller r s , with the situation 
of a heterojunction with a metal gate which screens the inter-electron Coulomb interaction, 
changing it to V e s(r — > oo) ~ 1/r 3 and V e s(r — > 0) ~ 1/r. For the gated heterojunction 
case clearly the Laughlin liquid will be preferred because the effective interaction remains 
Coulombic for small r. Thus, careful experiments in wide wells contrasted with those in 
gated heterojunctions would go a long way in establishing the liquid-solid phase boundary 
in strong-field 2D systems. 

The authors wish to thank P. I. Tamborenea and Song He for helpful discussions, and 
P. I. Tamborenea for providing us with his data. This work was supported by the US-ARO, 
the US-ONR, and NSF. 
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FIGURES 

FIG. 1. The effective interaction (solid line), Coulomb interaction (dash-dot line), and model 
interaction (dashed line) for a square well of width L. The model interaction parameter A = 0.2L. 
The lowest subband cosine wavefunction and gaussian fit are shown in the inset. 

FIG. 2. The variational liquid energy (dashed line) and solid energy (solid line) shown after 
subtracting %uj c /2 and the Madelung energy for various fixed values of r s . 

FIG. 3. The variational liquid energy (dashed line) and solid energy (triangles) shown after 
subtracting huj c /2 and the Madelung energy for various fixed values of A. 
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